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ABSTRACT 

We explore the pulsar slot gap (SG) electrodynamics up to very high altitudes, 
where for most relatively rapidly rotating pulsars both the standard small-angle 
approximation and the assumption that the magnetic field lines are ideal stream 
lines break down. We address the importance of the electrodynamic conditions at 
the SG boundaries and the occurrence of a steady-state drift of charged particles 
across the SG field lines at very high altitudes. These boundary conditions and 
the deviation of particle trajectories from stream lines determine the asymptotic 
behavior of the scalar potential at all radii from the polar cap (PC) to near 
the light cylinder. As a result, we demonstrate that the steady-state accelerating 
electric field, E\\ , must approach a small and constant value at high altitude above 
the PC. This is capable of maintaining electrons moving with high Lorentz 
factors (~ a few x 10^) and emitting curvature 7-ray photons up to nearly 
the light cylinder. By numerical simulations, we show that primary electrons 
accelerating from the PC surface to high altitude in the SG along the outer 
edge of the open field region will form caustic emission patterns on the trailing 
dipole field lines. Acceleration and emission in such an extended SG may form the 
physical basis of a model that can successfully reproduce some pulsar high-energy 
light curves. 

Subject headings: acceleration of particles — gamma rays: theory — pulsars: 
general — radiation mechanisms: nonthermal — stars: neutron 



1. INTRODUCTION 

There is no doubt that pulsars are accelerating particles up to relativistic energies in their 
magnetospheres, and that these particles are primarily responsible for the pulsar radio- to 
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high-energy non-thermal emission. It is also believed that the energetics of this acceleration, 
as well as the main physical processes involved in production of high-energy photons, are more 
or less understood. However, the ambiguity in interpretation of pulsar timing observations 
in terms of emission site mapping in a pulsar magnetosphere makes it difficult to answer the 
basic question of where the pulsar high-energy emission originates. In their recent attempt to 
explain the observed high-energy light curves of pulsars, Dyks & Rudak (2003) concentrated 
on a purely geometrical model by postulating that the emission is produced in a relatively 
narrow region along the last open magnetic field lines of a pulsar magnetosphere. The 
interesting result of their study is the occurrence of caustic emission zones (Morini 1983), 
i.e. the phase shifts of radiation emitted at radii between ~ 0.1 — 0.7 times the light cylinder 
radius, parallel to field lines on the trailing edge of the polar cap (PC), are cancelled by 
phase shifts due to rclativistic effects of aberration and time- of- flight. Radiation emitted 
over a large range of altitudes thus arrives in phase, forming two narrow peaks in the light 
curves, very similar to those of known 7-ray pulsars (e.g. Thompson 2001). 

In our previous paper (Muslimov & Harding 2003 [MH03]) we began discussing the 
regime of acceleration of particles and production of high-energy emission within the pulsar 
slot gap (SG), a narrow region on the boundary of the open field lines, where the electric field 
drops to zero. The SG is a pair-free region of slower acceleration, in which the parallel electric 
field is unscreened. Pair cascades develop along the inner edge of the SG at several stellar 
radii above the NS surface. Even though the SG regime in pulsars was originally introduced 
in the electrodynamic model of Arons & Scharlemann (1979), it was not considered a viable 
high-energy emission region (see e.g. Arons 1996). The revised version of the SG regime 
proposed by MH03 incorporates the effect of rclativistic frame dragging (Muslimov & Tsygan 
1992 [MT92]) and, more importantly, the effect of SG boundaries on the strength of the 
accelerating electric field within the SG. MH03 demonstrated that the primary electrons 
tend to accelerate up to higher altitudes before pair production begins, and pair cascades 
continue along the inner boundary of the SG until the magnetic field becomes too low. The 
resulting radiation from the pair cascades forms a wide, hollow cone of high-energy radiation 
due to the flaring of fleld fines. Adhering to the small-angle approximation, MH03 restricted 
their study of the SG regime to altitudes less than four-flve stellar radii. However, since the 
parallel electric fleld in the SG is not screened on fleld lines close to the open-fleld boundary, 
acceleration may continue to much higher altitudes. Particle acceleration and radiation in 
such an extended SG may therefore provide a physical basis for the two-pole caustic model 
of Dyks & Rudak (2003). 

Formation of a SG requires the production of enough pair multiplicity to screen the 
parallel electric fleld above the pair formation front. We have found from our previous studies 
(Harding & Mushmov 2001 [HMOl], 2002 [HM02]) that the youngest and most energetic 
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pulsars can produce pairs from curvature radiation (CR) of primary electrons, which are 
numerous enough to screen the electric field. Older, less energetic pulsars, those below 
the CR pair death line, can produce only pairs from inverse Compton radiation of primary 
electrons scattering thermal X-rays from the NS surface. The inverse Compton pairs are not 
numerous enough to completely screen the parallel electric field. A necessary condition for 
formation of a SG is thus the ability to produce pairs from CR, and the expression for the 
CR death fine (given by Eqn [52] of HM02) defines the boundary in the P-P diagram of 
pulsars capable of having SGs. Such pulsars include the Crab, Vela, Geminga and most of 
the 7-ray pulsars detected by EGRET, but not the majority of miUisecond pulsars. 

The extension of the regime of SG acceleration to much higher altitudes is the main 
subject of the present paper. In the Sections below we outline our approach to constructing 
an appropriate steady-state physical solution that can be used up to very high altitudes in 
the SG. We also discuss the immediate consequences of our proposed extended SG solution: 
acceleration of particles (electrons, positrons) and high-energy emission up to nearly light- 
cylinder radius, and the possibility of occurrence of high-altitude caustic emission on trailing 
field lines. 

The paper is organized as follows. In §2 we discuss the electrodynamics within the SG 

regions of pulsars. We address the physical constraints on the scalar potential (§2.1) and 
equipotcntiality of SG boundaries and derivation of effective Poisson's equation (§2.2) in 
the outermost section of SG. In §3 wc present the electrodynamic sohition within the SG at 
very high altitudes. In §3.1 we illustrate the possibility of extended acceleration within SG 
in the regime where the acceleration is balanced by the curvature-radiation reaction. Our 
numerical calculations are discussed in §3.2. Finally, in §4 we discuss our main results and 
draw our principal conclusions. 



2. Steady State SG Electrodynamics 

In the frame of reference rigidly corotating with a neutron star (NS), where the magnetic 
field is stationary (and having pure dipolar geometry), the general relativistic Maxwell's 
equations yield (see MT92) 

E- — (w-u) X B = --V$, (1) 
ac a 

where E and B are the electric and magnetic fields defined in Zero- Angular-Momentum- 
Observer (ZAMO) frame of reference (see Macdonald & Thorne 1982), u is the rotational 
velocity, and w is the differential velocity of rotation of inertial frame of reference, $ is the 
scalar potential, and the so-called geneneral-relativistic 'lapse function', a, is defined below. 
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right after expression (18). Taking the divergence of eq. (1) and making use of Maxwell 
equation 

V • E = 47rp, (2) 
we get the Poisson's equation for the scalar potential $ 

-V^) =47r(p,,-p), (3) 



a 



where 



— (u - w) X B 

a 



is the general relativistic expression for the Goldreich-Juhan (GJ) charge density (cf. Gol- 
dreich & Julian, 1969), and p is the actual charge density of electrons determined by their 
relativistic flow along the magnetic field lines and which is fixed by the condition E\\ = at 
the stellar surface (see MT92 for details). Thus our electro dynamic description of charges 
streaming along the open field lines will imply the space-charge-limited flow approximation 
(at least near the stellar surface, within the radial distance of less than a few stellar radii). 

Note that the l.h.s. of eq. (1) can be treated as the effective electric field in the frame 
of reference rigidly corotating with the NS, 

E' = E- — (w-u) X B. (5) 
ac 

The condition of absence of any electric field in the regions of the magnetosphere with closed 
field fines rigidly corotating with the NS, E' = 0, necessarily implies that these regions should 
be filled with charges of density p — p^j (as a trivial solution of eq. [3]). 

Our previously derived solutions (see MT92, Muslimov & Harding 1997 [MH97]) for the 
case IpI < jp^jl were limited by a small-angle approximation and therefore cannot be justi- 
fiably used beyond the radial distances of about ~ 3-4 stellar radii above the PC surface 
(for a typical pulsar spin period). Recently, MHOS discussed the SG solution which is also 
formally limited to a small-angle approximation. This means that, because of the curving of 
the SG toward the magnetic equator, the solution derived in MHOS cannot be used for high 
altitudes, typically exceeding a few (or several, at most) stellar radii above the surface. Here 
we discuss the regime of steady-state acceleration within the SG extending up to very high 
altitudes, nearly approaching the light-cylinder. By addressing the basic physical conditions 
that are required for the occurrence of this regime, we observe that the standard concept of 
GJ charge density becomes inapplicable in the outermost section of the SG, where the effec- 
tive GJ charge density gets significantly constrained by the requirement of equipotentiality 
of SG boundaries and by the effect of cross-field motion of charges. 
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We shall explore the SG solution in the outermost part of NS magnetosphere (but still 
within the light-cylinder) satisfying the same boundary conditions as those used in all our 
previous studies: equipotentiality of the SG surfaces, and zero-electric field condition at the 
PG surface. We propose the following ansatz for constructing the general solution extending 
from the innermost section up through the outermost section of the SG. By considering 
the property (see Section 2.1 below) of the scalar potential $ at large (up to the light- 
cyhnder radius) distances found by Mestel et al. (1985) together with the abovementioned 
boundary conditions (equipotentiality of SG surfaces and zero-electric field condition at the 
PG surface), we can unambiguosly constrain and determine the outermost solution for the 
potential Then, by matching the outermost solution with the known innermost solution 
near the NS surface (presented in MHOS) wc can construct the approximate general solution 
applicable to both the innermost and outermost sections of the SG. 

Note that any physically meaningful electrodynamic solution in the outermost part of 
the magnetosphere should take into account the effect of particle drift across the field lines 
or deviating of particle trajectories from magnetic field lines. This effect should unavoid- 
ably constrain the scalar potential simply because in this region the field lines cannot be 
treated as characteristics or as stream lines for the flux of relativistically moving electrons. 
In a steady-state situation, it is reasonable to expect that the scalar potential $ is a mono- 
tonically increasing (or decreasing, as in the case of acceleration of positive charges) and then 
saturating function of radial distance so that the outermost solution gradually matches the 
innermost one. In this study we demonstrate that in a steady-state situation, the constraint 
on $ in the outermost region of the NS magnetosphere (but well within the fight cylinder) 
along with the condition of equipotentiality of the SG surface allows us to derive an appro- 
priate electrostatic solution. This solution implies initial (in the innermost section of a SG) 
boosting of electron acceleration over characteristic lenghtscale of ~ 1 — 2 stellar radii and 
subsequent extremely slow post-boost acceleration over lengthscale ~ light- cylinder radius. 

2.1. Constraint on Potential $ in the Outermost Section of SG 

The equation of motion of an electron of mass m and charge -e can be written as (see 
eq. [2.24] in Mestel et al. 1985) 




(6) 



where 7 = (1 — /c^) v is the electron velocity. In Mestel et al.'s notations. 



E = -V0, 



(7) 
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and 

I X f B - —V X m) = V(/)*, (8) 
c \ e / 

so that equation of motion, (6), translates into (see eq. [2.26] in Mestel et al. 1985) 

-e(f) + -fmc^ ^ -e(f)*{S), (9) 

where S is the stream function. 

Note that the scalar potential $ defined by eq. (1) is a general- relativistic counterpart 
of the so-called 'non-corotational' potential ip (see e.g. eq. [2.5] in Mestel et al. 1985) 
introduced by Endean (1974), Mestel (1973), and Westfold (1981) and which can be defined 
via equation 

VV^ = V(/) - X B. (10) 

For relatively low altitudes $ ~ 7(mc^/e). For high altitudes (but still within the light- 
cylinder), as was first demonstrated by Mestel et al. (1985), the change in the angular 
momentum ^ymftr^ (where f is the radial cylindrical polar coordinate) of a streaming particle 
occurs only through the toroidal component of the magnetic force, requiring departure from 
strict flow along the field lines. The combination of energy and angular momentum integrals 
gives (in cylindrical polar coordinates in the axisymmetric case) 

^ = ^,{i-^)-ns), (11) 

where T{S) is some function which is constant on stream lines and which is set at the stellar 
surface. Here the term cx fl^ is the so-called 'centrifugal-slingshot' term (see Mestel et al. 
1985) arising from the fiow of electrons across the field lines. 

It is important to point out that within the domain of the SG the actual deviation from 
strict flow along (poloidal) B is of order 

5^^mcy{eBAl^^), (12) 

where AZg^ ~ ^o-R^{sg\A/* characteristic latitudinal SG thickness at dimensionless 

radial distance rj^ [r] = r/R); 9q is the PC half-angle, = [^-R//(l)c]^^^; and Ai^g^ is 
the latitudinal SG thickness in units of ^{^ — 9 / 9q \s the dimensionless colatitude of a PC 
field; see Section 3 and also MH03 for details). Thus, the characteristic dimensionless radial 
distance at which the magnitude of the deviation from strict fiow along (poloidal) B reaches 
5, can be estimated as 

~0.6,..(A4,5)^/^^(|^)''', (13) 
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where rjic — c/QR is the dimensionless radius of the hght-cyhnder; Q = 27r/P is the angular 

velocity of NS rotation; R and P are the NS radius and spin period, respectively; B12 — 
Bq/IO^'^ G, Bq is the surface value of NS magnetic field strength; Rq — R/10^ cm, Pq.i — 
P/0.1 s, and 77 = 7/10^. 

For the parameters of the Crab pulsar (here we adopt the values: B12 — S, Rq — 1.6, and 
Po.i = 0.33) and assuming 5 ~ 0.05 — 0.1, from formula (13) we get 

,. = (0.3-0.4),. (^1)"'. (14) 

SO that for the estimated value of A^g^ ~ 0.05 (see MH03) and for 77 ~ 3 condition (13) is 
satisfied at r]^ ~ (0.5 — 0.7) rjic- 

Formula (14) means that within the SG the effect of transfield motion becomes important 
already at rj<r]ic (for the Crab-like pulsars) and should be taken into account. This is 
significantly different from Mestel et al.'s model where a similar situation would occur well 
beyond the light cylinder (and where A/g^ in formula [12] should be replaced by the light 
cylinder radius, Ric = c/fl), because the lengthscale was the entire open field region rather 
than the narrow SG. 

The fundamental consequence of eq. (11) is that in the region where the centrifugal- 
slingshot effect becomes important (e.g. at < r/ic foi' the Crab-like pulsars) the electrons 
crossing the field lines begin picking up energy from the corotational part of the potential 
(from potential 0, in Mestel et al.'s notation), so that in the regime of mostly transverse 
flow the change in potential (change in T{S)) across the field lines caused by rotation tends 
to balance the change in the corotational part of the potential. Suppose that the electrons 
are flowing with relativistic velocities along the magnetic field lines and entering the region 
where they are getting 'decoupled' from the magnetic field lines. Apparently, the solution 
of MH03 for p and therefore for $ will not be warranted in this region and especially in 
the region with predominantly transfield flow. However, we may justifiably assume that at 
the onset of the transfield flow regime, where the relativistic flow of electrons is still mostly 
along the magnetic field lines, the condition (see e.g. Mestel 1995, 1999) 

E'x ^E± + X B 0, (15) 

c 

'turns on', with E\\ ^ E± (for most acceleration scenarios we discuss in this paper E\\ is 
balanced by the CR reaction force). Further out in the magnetosphere the condition (15) 
may transform into the perfect MHD condition (see also Contopoulos, Kazanas & Fendt, 
1999), E -|- (v/c) X B = 0, which is a good approximation as long as E\\ -C E±. 
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We suggest that the condition (15) together with the equipotentiahty of the SG surfaces 
may in fact determine the behavior of potential $ through the outermost section of the SG. 
It is important that at large radial distances E'_|_ as given by eq. (5) is dominated by the term 
(O X r/c) X B, which tends to produce significant electric potential drop across the SG, and 
as a result, would induce enormous surface charge on the SG boundaries. Also, it must be 
pointed out, that this term would tend to induce a strong component of i^y with the polarity 
that may change with altitude and become opposite to that of the main component of E\\, 
produced by a small imbalance between the GJ and actual charge densities in the innermost 
region, at altitudes within ~ 1 — 2 stellar radii. The occurrence of a strong component of 
£■11 with reversed polarity would unavoidably disrupt the continuous flow of electrons along 
the field lines and result in an essentially non-stationary regime of particle flow. However, 
in the steady-state situation we consider in this paper, the occurrence of cross-field motion 
of electrons would effectively screen out the excessive GJ space charge and short out the SG 
boundaries, thus maintaining them as equipotential. The latter means that condition (15) 
would be roughly satisfied in and beyond this region. Before we discuss how condition (15) 
can be explicitly incorporated into our electrostatic solution, let us discuss the consistency 
of this condition with the assumption of equipotentiahty of SG boundaries all the way from 
the PG surface up to the very high altitudes (say, up to ~ 0.1 — 0.5 of the hght-cylinder 
radius) . 



2.2. Equipotentiality of SG Boundaries and Effective Poisson's Equation in 

the Outermost Section of SG 



Let us consider the cross-sectional area of a magnetic flux tube (of dipole field) emanat- 
ing from the PG at radial distance rj {— r/R) 



'5(^)='5(l)ffi|,^ 

where /(r/) is the general-relativistic correction factor, defined as (see e.g. MT92) 



/(^) = -3 ( 



In 1 - 



+ 



1 + 



2r] 



(16) 



(17) 



where r) — r/R is the dimensionless radial coordinate, e — Tg/R, and is the gravitational 
radius of the NS. One can use approximate formula /(?]) ~ l-|-0.75a;-|-0.6a;^ (where x — s/rj). 
For a canonical NS of 1.4 solar mass and 10 km radius {e — 0.4) /(I) 1.4. 

In this paper, as in our previous studies, we use the magnetic spherical polar coordinates 
{rj,9,(f)pc)- We also denote by x the pulsar obhquity (angle between the NS rotation axis 
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and magnetic moment). We will refer a 'normal polarity' pulsar as one having 0° < x < 90° 
(north magnetic pole near north astrographic pole: f2 • m > 0, where m is the NS magnetic 
dipole moment), and a 'reversed polarity' pulsar as one having 90° < x < 180° (north 
magnetic pole near south astrographic pole: O • m < 0). 

We shall now introduce the flux of charges streaming with relativistic velocity through 5(77), 

F = a{i^)cS{i^)p{i^), (18) 

where p is the local charge density, and a — {1 — e/rj)^/'^ is the lapse function. 

In a steady-state regime, well within the light cylinder, the flux F should be constant along 
the individual magnetic flux tube as a consequence of charge continuity equation which 
implies that ap oc B cc r]~^. In this case F is a function of ^ only (and not 77). [The variable 
^ is equivalent to the Stokes stream function S (see eq. [9]) used by Mestel et al.] However, 
at very high altitudes, for any stream line within the SG we can write that 

B = (19) 



where f3 = \J\ — Srj/Arjic is a factor that takes into account the change in the geometry of 
the flux tube as we move from the magnetic pole to the equator. 

The explicit expressions (see e.g. MH97, for the derivation of these expressions) for the 
GJ charge density, p^j, and actual charge density, p, may be written as 

f f 1 1 

Pgj = ~^077TT — §{K(0 + «i(^>0]cosx + [&o(0 + ^i(^>O]sinxcos0pc}, (20) 

and 

f(v) 1 

P ~ -^^0^^^/3(77) {[ao(0 + 0] cosx + MO + hih 01 sinxcos0pe} , (21) 

respectively (see also MT92 for the exact expression for p valid for arbitrary small altitudes 

and which is consistent with the radial profiles of potential $ and E^y). Here po — flBo/2Trc, 

and ^ g 

ao = 1, ai = -4 - ::H{7]) sin^ 9, (22) 

77 z 

60 = 0, bi = ^H{r]) sin 9 cos 9, (23) 

-r-r, S ^ >^ 1 3 £ K \ _ 
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Here n is the general-relativistic parameter characterizing the magnitude of the frame- 
dragging effect near the stellar surface measured in stellar rotation velocity, f2, and /3 is 
defined right after eq. (19). These expressions are formally derived for arbitrarily large 
distances but still within the light cylinder. Well within the light cylinder, at <^ r]ic 
{flic ^ (^o^)-i ^-iid in a small-angle approximation, 

ai--^, ^^-^lWri)p^y (25) 

which correspond to the expressions for p and p^j we used in our previous papers for the 
situations where a small-angle approximation was in fact more than satisfactory. One can also 
use the approximate expression, H{ri) ^ 1 — 0.25x — 0.16x^ — 0.5(«;/£^)x^(l — 0.25x — 0.21x^), 
where x — e/rj. For a canonical NS [see also eq. (17)] H{1) 0.8. 

Now, the flux of effective GJ charges can be written as a sum of constant and varying 
with altitude components, 

^G.=^G.,0+^G..l, (26) 

where 

-^Gj,i = -^o(ajcosx + &isinxcos0pc), i = 0, 1 , (27) 

and Fq — —pocS{l) is a constant factor depending on pulsar bulk parameters. Note that, 
according to the GJ reasoning (Goldreich & Julian, 1969), the corotating region of the pulsar 
magnetosphcrc with closed magnetic field lines should be filled with charges of local density 
Pqj. This condition guarantees that any electric field which could be possibly generated 
in this region should be completely screened out. We do not intend to speculate on the 
dynamics of formation of the non-vacuum pulsar magnctosphere (see e.g. Krause-Polstorff 
& Michel, 1985; Arons & Spitkovsky, 2002), but we assume that the filling in of the closed 
field fines of the initially charge-starved magnctosphere with charges is most likely to occur 
along the field lines. Even more, since F^j ^ is a function of rj, one may expect that this 
process develops in a non-stationary manner. 

Let us now consider the flux AF, corresponding to the charge imbalance between the 
GJ and actual local space-charge density, Ap = , — p, 

AF = F^^-F^ acS{rj){p^, - p) ^ AF^{r],0 - /3(7;)AF,(1, 0, (28) 

where 

A^*(?7,0 = ^oK(?7,0 cos X + &1 (77,0 sin X cos 0pc] (29) 
Using the expficit expressions for ai and bi (see eqs [22], [23]) we can write 



AF a; Fq 



cosx + 
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^^(1)^0 I - /3 1 sinxcos^pe } , (30) 



2 

where /? is defined right after eq. (19). Here we should reiterate that the formal usage 
of a general expression for Ap = p^, — p (= AF/acS[ri]) to solve the Poisson's equation 
for arbitrarily large altitudes leads to the inconsistent and even erroneous result, mostly 
because of the physical reasons discussed in the end of previous Section. Namely, such 
a solution would imply the building up of the effective flux, AF, exceeding the GJ flux 
Fqj(77 = 1), fixed at the stellar surface, and sign reversal of the accelerating electric field 
at high altitudes. Also, we would like briefiy comment on the space-charge-limited fiow 
approximation at high altitudes in pulsars. Generally, the space-charge limitation occurs 
when the ejected charges reduce the accelerating potential drop boosting the initial particle 
energy. In pulsars, the flux of electrons ejected from the PC surface is limited by the value 
of GJ space charge at the bottom of the PC. In this case F remains constant along the 
magnetic stream lines and is determined by F^^ at rj — 1. Above the PC surface the space 
charge of ejected electrons reduces the "vacuum" potential drop by limiting it to the value 
determined by a small imbalance between the GJ charge density and actual charge density 
of electrons. This approximation is perfectly valid within a few stellar radii above the PC 
surface of most pulsars, when |AF(r;)| < |F(-,j(l)|. However, at high enough altitudes, where 
|AF(?7)| > |Fqj(1)| the situation is akin to the acceleration of test particles in a vacuum-like 
potential drop. Let us examine this in more detail for the SG by using the above expressions 
for AF and then formulate the derivation of an approximate but physically meaningful 
solution. 

For low altitudes, 77~l-|-^(^-Cl),we can write 

AF 3Fo^ (^K cos X + ^H{l)eo sin x cos 4>p^ , (31) 



k) cos X + 7;H{l)9o sin x cos 0pc 



(32) 



and 

3 
2 

so that I AF| <^ l-^GjI' assuming that F^^ ^ 0. This means that at low altitudes the imbalance 
between F and (the flux of flctitious charges) that gives rise to the accelerating electric 
fleld in that region is much smaller than the local GJ flux, and therefore there will be no 
disruption of the steady-state regime of particle flow within the SG. 

At large distances, 1 <S 77 <S r\ic-, 

|AF|^|AF,(r;,OI«^|i^Gj.il, (33) 
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where I-Fqj 1(77 » 1)| ^ I-^gj(-^ = ^ ~ 1 ^ 1)L ^-iid -Fqj(2; <^ 1) is a function of ^ only (i.e. is 
nearly a constant along the field hnes) . 

Thus, for low altitudes, up to approximately one-two stellar radii above the PC surface, 

^GJ~^GJ,0, (34) 

and 

< \FJ. (35) 

Expressions (34), (35) imply that the SG boundaries can easily acquire the necessary surface 
charge by means of redistribution of space charges along the field lines in the vicinity of the 
SG boundaries, thus enabhng the fulfillment of both the equipotentiality of SG boundaries 
and continuity of across them. Note that, the component and therefore the magnitude 
of induced surface charge is mostly determined by -Fqj ^ (or by PqjJ- Since in this region 
I-P'gjiI — l-^GTol' required surface charge can be easily built up by establishing a weak 
current along the boundary field lines, without violating the GJ condition that in a steady- 
state regime the fiux of charges from (to) the stellar surface should be limited by F^,,,,. We 
must also note that the presence of a weak current along the SG boundary determined by 
flux j is perfectly compatible with equipotentiality of the boundary. For example, this 
can be achieved by establishing a slightly non-homogeneous distribution of charges in a tiny 
skin layer along the boundary, at $ = 0. 

As we move up to higher altitudes along the SG boundaries where 

\FaJ>\FoJ, (36) 

and where E'^ is mostly determined by the term ~ (fi x r/c) x B, the situation changes 
dramatically. In this case the equipotentiality of the SG boundaries becomes fundamentally 
incompatible with a steady state regime. In other words, it is very unlikely that the SG 
boundaries can be steadily maintained in dynamic equilibrium in the presence of a strong 
E'j^ component. Rather, it is this region where the centrifugal-shngshot effect makes the 
electrons/positrons 'slip' from the magnetic field lines, and therefore effectively prevents 
charges from building up the otherwise required surplus surface charge on SG boundaries. 
Thus, in a steady-state situation the SG boundaries can be maintained as equipotcntial, if 
the approximate condition E'_l = (or its classical counterpart [15], as it will be referred 
to in the rest of the paper) is achieved throughout the SG outermost section. This means 
that in the outermost section of SG the value of AF(?], ^) (see [30]) cannot grow (because of 
cross-field motion of charges at very high altitudes that effectively destroys the excessive GJ 
space charge in this region), and it is likely to nearly saturate at 77 ~ 77c remaining constant 
along the stream hnes. In this case the scalar potential $ is described by the following 
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Poisson's equation 



a 



V$ = 47rApeff, 



(37) 



where ApeS is now determined by AF{r)c, (see formulae [20], [21], [28] and [30]), and reads 

m 1 



'eff 



-Po 



/(I) ar)^ 



{[1 -/3 + ai(?7c,0 -«i(l,0/5] COSX + [6i(?7c,0 - &i(l, 0/^] sinx cos 0pc} 



-Po 



/MJ_ 

/(I) ar^3 



1 



M/?-;g)+l-/5 



COSX+ 2^(1)^0 



HjVc) hem 

H{1) V /(^c) 



- (5 sinxcos^pc >(38) 



where (3 is defined after eq. (19). The parameter rj^ entering the r.h.s of eq. (38) and defining 
the altitude of saturation is to be determined through matching the solution at low altitudes 
of MH03 (see Section 3.2). 



3. Extended-SG Electrodynamic Solution 

The general formula for the dipolar magnetic field line within the SG can be written as 
(see MH03) 



'^y|^^o,sc(lT^Ae3oe*), (39) 

where 9^ ^^ = 9o{l - A^gQ/2) ^ = [^^/c/(l)]^''^ is the polar angle of the SG central line. 
This central line separates the SG into the innermost half-space and the outermost half- 
space where the dimensionless colatitude of field lines C,* varies from 1 to and from to 1, 
respectively. Here A^^q is the latitudinal gap thickness in units of dimensionless colatitude 
i a = 6'/6'o). Note that in a small-angle approximation one may justifiably assume (as it 
was done in HM03) that sin ^ fa ^ and 6'q?7 <^ 1. 

The above formula for the magnetic field lines is based on a canonical definition of the 
light-cylinder radius, Ric = c/Q. Strictly speaking, this definition is valid only for a nearly 
aligned rotator with poloidal dipolar magnetic field. For a NS with arbitrary obliquity, X) 
the footpoints of last open field fines slightly deviate from those in the case of an ahgned 
rotator, and the magnetic colatitude of the last open field line becomes a function of x 
and 0pc. In addition, the rotational distortion of a pure dipolar poloidal magnetic field 
near the hght-cylinder also affects the shape of the surface formed by the last open field 
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lines and location of their footpoints. The non-circularity of the PC boundary caused by 
pulsar obliquity and rotation were studied by many authors (see e.g. Cheng, Ruderman & 
Zhang, 2000 and Dyks, Harding & Rudak, 2004, for most recent assessment of this issue; and 
references therein). It is interesting, that for x<(60° — 70°) the assumption that the light- 
cylinder radius is simply defined as c/r2 seems to be very satisfactory in most theoretical 
studies on pulsar magnetospheres. For example, for x ~ (60° — 70°) the effects of pulsar 
obliquity and magnetosphere rotation discussed in Cheng, Ruderman & Zhang (2000), may 
distort the circularity of an "ideal" (very small obliquity and no rotational distortion) PC 
boundary by less than 20-25%. 

We shall reasonably assume that A^g^ -C 1 (thin SG approximation). For \ <^rf <^ rjf^ 
(where rjic — Ric/R), the operator in the l.h.s. of Poisson's equation (see eqn. [37]) 
reduces to 



sin e 09 \ 09 J sm^9 0(t). 

SI fM^^ ' (40) 



Changing the variable 9 to (via relationship [39]), we get 



2 

pc 



where z/g^, = (1/4) A.^^^ (sec MH03 for more details). Thus, by using expression (41) the 
Poisson's equation (37) can be rewritten as 

1 0'^ f(l) 1 0'^ \ 

— — ± ; v9l^— + -— U = -2^o9l^^{Acosx + Bsmxcos(Pp,), (42) 



* ^ ^ /43N 

- 1-vfiiKJfiv) ^ ' 



where $o = {^R/c)BqR, and 



Here A and S, according to eq. (38), are given by 



.4 = .(/3-y+l-„ and B=5W)(^^^-/3|, (44) 

respectively. In a small-angle approximation (see MH03) one would neglect the second term 
in the l.h.s. of eq. (42) and replace v*^ by v^^. 
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We shall search for general solution of eq. (42) in the form 

$ = $1 cos X + $2 sin X cos (f)pc- 
The general solutions for $i and $2 read 
2/(^) 



$1 = ±- 



and 



cosh(y^{* 



$2 = '^^e^.^B 



1 - 



cosh(^/i^) 



respectively. 

In expression (46), 



/(I) 



SG 



7(^)1 -^/(l)^o%c//(^)' 



When s < 1 (or 1 - 77^2 ^ 1 - r]/r]ic > ^/i^), 

^ ^ SG 0,SG 

Thus, for 1 — rj/rjic^ ^li® above solutions for $1 and $2 yield 



$(77, = $o<3G ^s*G^(l - cosx + 2B 



1 - 



cosh(yz^^^) 



cosh(^/I^) 

For 1 -C 77 < eryjc (where e ~ 0.1 — 0.5), eq. (50) can be rewritten as 



sin X cos 0^ 



pc 



V c J /(I) 



H{i) V'm 



p 



1 -C*)cosx+ 
cosh(^/z^{,.)" 



cosh(yi^) 



sm X cos 



(45) 



(46) 



(47) 



(48) 



(49) 



(50) 



(51) 



In a thin-SG approximation {^/i^ <^ 1), the second term in eq. (51) can be simplified to 
yield 



- 



«:(/3-^)+l-/3 



1 H ) cos x+ 

Vic. 



H{1) V /(^c) 



Vc^TZ^ - I sinxcos0pc (1 - C*)- 



(52) 
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This solution illustrates also that for each magnetic stream hne (each ^) within the SG the 
scalar potential $ nearly saturates at high altitudes (but well within the light cylinder). The 
remarkable property of the SG solution is that it incorporates the effect of the SG bending 
(/5 and 1 + rj/rjic terms in formula [52]). The first (/j-) effect is owing to the dependence of 
p on at high altitudes discussed in Section 2.2. The second effect is due to the fact that, 
as the altitude increases, the divergence of the electric field within the SG tends to decrease 
because of the straightening of the SG (see the first term in eq. [42]). However, since the 
charge per unit length of the SG remains constant (see r.h.s. of eq.[42]), this geometrical 
effect is compensated by slight enhancement of potential $ at higher altitudes. Finally, 
solution (52) reduces to the low-altitude solution derived in MH03 if we neglect in formula 
(52) the term rj/rjic and replace r]c by rj. Note that in this case $ = — (2/c)AF(l — ^^), i.e. 
the scalar potential $ is simply proportional to the flux of flctitious charge, AF. 

Before we discuss how to determine parameter r]c entering solution (52) we shall consider 
more closely the regions with the so-called "favorably" and "unfavorably" curved field lines 
introduced by Arons & Scharlemann (1979). The region with favorably curved field lines 
simply corresponds to that with cos 0pc > (electron acceleration, for normal polarity), and 
the region with unfavorably curved field lines corresponds to that with cos < (positive 
charge acceleration). Note that for the reversed polarity pulsar (see definition below formula 
[17]) the sign of accelerating particles reverses. Since in the electrodynamic model of Arons 
& Scharlemann (1979) the accelerating electric field is solely determined by the component 
oc sin X cos (f)pci the electrons or positive charges (e.g. positrons) may accelerate in their model 
only along favorably or unfavorably curved field lines, respectively. In this sense, in our model 
the presence of the component oc cos x (which is also independent of 0pc) tends to make all 
field lines 'favorably curved' in that the same sign of charge is accelerated outward over the 
entire PC. However, for very large obliquities and/or very short spin periods the favorably 
curved field lines may again become unfavorable for electron acceleration, because of the 
negative contribution from the term oc sin^cos^pc (for the region where cos^pc < 0). Now 
consider the situation where in the formula for $ the component which is oc cos x dominates. 
In this case the value of r]c can be estimated from matching the asymptotic value of E\\ with 
the value of E\\ generated at small altitudes, which translates into Sn/r]^ {3 + 5K,)/Sriic. For 
the parameters of the Crab pulsar and assuming very small obliquity this condition would 
result in rjc ~ 3. For x — 60° and cos (ppc < 0, we come up with the estimate rjc ~ 1.2 — 1.3. 

Using formula (52) and E\]^ — — (1/q;)V||$ = —{l/R)d^/dr) we can get, for rj ~ rjc, 



3 /nn 




3 



Bo 



El 



,high 



/(I) 



SG 



V 

— COS x+ 

Vlcj 




(53) 
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According to formula (53), 



^ii,high| ~ ::§:VsG < ^(^ ~ vic) ~ \E±{r}i^)l (54) 



Vic 

which is the condition accompanying the regime given by Eqn (15). 



3.1. Extended Acceleration Within the SG in Curvature- Radiation Reaction 

Limit 

Let us estimate the characteristic electron Lorentz factor, 7, in the regime of extended 
acceleration within the SG in the situation where the gain in electron kinetic energy is almost 
exactly compensated by the CR reaction, i.e. where the condition 

2 

e^ll,high ~ ^—7^ (55) 

is satisfied. Here Pc-^ {^l?)) R {'nrii^^/'^ is the characteristic curvature radius of the field hnes 
in the SG. 

By using formula (53) in the above equation, we arrive at 



7 ~ 3 • 10^ 



o3 „ ni/4 
(2.5 + 0.6«o.i5)Si2^— i^sG(l - cosxl 

^0.1 Vic 



(56) 



where Ko.15 = fi:/0.15 ^5 = I/IO^^ g-cm^, I is the NS moment of inertia. Here, for 

the sake of simplicity, we assume that the term oc cos x dominates. 

Using parameters of the Grab pulsar, B12 — 8, Pq.i — 0.33, Rq — 1.6, and for ~ 0.05 and 
77 ~ O.Brjic we get 

7 ~ 4 • 10^[(2.5 + 0.6«o.i5)(l - cosxl]'/'. (57) 

This estimate implies that the energetics of the extended SG region is about the same as that 
of the region near the PG surface, even though the physical conditions for pair production 
and generation of high-energy emission in the outermost part of the SG are quite different 
(we shall address this issue in the next Section). 



3.2. Numerical Results 



In this section we will present some results of numerical simulations of the extended 
SG solution derived earlier in Section 3 for a particular case. The results shown in this 



-18- 



Section are meant to illustrate the most important features of the proposed acceleration 
model. A more extensive exploration of pulse profile properties and spectra, and fitting to 
observed high-energy pulsars, will appear in another paper. We wish to have an approximate 
expression for the accelerating E\\ in the SG from the PC surface to near the light cylinder. 
As discussed earlier in Section 3, we are not able to derive a single solution valid at all 
altitudes, but were able to find expressions both within a few stellar radii of the PC surface 
and at large radii, when condition (15) dominates the solution. Our first step is thus to 
find an expression that best matches the solutions in the two regions, by finding the value 
of parameter rjc- We will then use this approximate expression to simulate acceleration and 
pair cascade radiation of primary electrons from the stellar surface. 



3.2.1. Matching Low Altitude and High Altitude Solutions for E\\ 

Eqn (12) of MH03 provides an expression for the scalar potential $ vahd in the SG at 
low altitudes (for radn rj<rjc). The low altitude accelerating Eh — — (1/q;)V||$ is thus 



£;ii,iow--3('^)'-^.3c(i-a 



K, 1 

— cosx + 7;^oH(l)5(r]) sinxcos0pc 



(58) 



where ^(77) = d\n[H6)/dr] (MH97) varies between ~ 0.5 and 1. One can write S{r)) = 
0.5 — d\nH/d\nx ^ 0.5 + 0.25a; + 0.39,t^ [where x = e/f]; see also eqs (17) and (24)], which, 
for a canonical NS, gives 5(1) ~ 0.7. As noted earlier, H{1) a; 0.8 for a canonical NS. 

At high altitude, rj ~ rjc, we have argued that the effective fiux, AF, will nearly saturate 
and that the high will be given by expression (53). In order to match these two solutions 
at high and low altitude, we write a general expression of the form 

E\\ ~ £;||,io^ exp [-{rj - l)/{r)c - 1)] + E^n.high (59) 

and determine 1]^ to give a smooth transition between -E||,high and -En^iow given by expres- 
sions (53) and (58), respectively. Formula (59) is valid (both for favorably and unfavorably 
curved field lines) when [k,cosx+ {^/'2)9oH{1)S{1) sin x cos (j)pc\ ^ or, more precisely, when 
AF{r])/Fo > (for r]<r]c). For a reversed polarity pulsar, formula (59) is also applicable for 
electron acceleration for nearly orthogonal rotator and for unfavorably curved field lines. 

For nearly aligned (0°<x<60°) or nearly anti-aligned (120° <%< 180°) cases, there will 
be continuous acceleration of electrons from the PC surface up to very high altitudes on both 
favorably (cos 4>pc > 0) and unfavorably (cos (f)pc < 0) curved field lines. For nearly orthog- 
onal rotators (80° <x < 100°), there will continuous acceleration of electrons/positrons from 
the PC surface up to very high altitudes on favorably curved field lines for normal/reversed 
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polarity, and continuous acceleration of positrons/electrons on unfavorably curved field lines 
for normal/reversed polarity. At intermediate inclination (60°<x~80° and 100° ~ 120°), 
there will be continuous acceleration of electrons on favorably/unfavorably curved field lines 
for normal/reversed polarity. However, on unfavorably/favorably curved field lines for nor- 
mal/reversed polarity, there may not be continuous acceleration of electrons from the surface, 
since the first term (oc cosx) in eq. (58) dominates only near the NS surface and the second 
term accelerates positrons. At higher altitudes, where eq. (53) applies, there will again be a 
regime favorable for electron acceleration. However, this situation is unstable since it implies 
a build-up of charge. The system will likely become self-limiting, possibly shutting down 
all charge flow along this held line. This situation may be favorable for the operation of an 
outer gap on such field lines (sec e.g. Cheng, Ho & Ruderman 1986). Thus, the extended SG 
may not provide steady-state particle acceleration at all inclinations. The value for rjc will 
vary with inclination angle and may approach e.g. close to the stellar surface for x ^ 60°. 

We have chosen for our numerical simulations an inclination angle x = 45° in order to 
illustrate the extended SG solution we discuss in this paper. For x — 45° we estimate that 
Vc = 1.4. 

3.2.2. Geometry of Acceleration and Emission in the Extended SG 

Using the approximate expression for E\i in eq. (59), we have simulated the acceleration 
and cascade radiation of a primary electron beginning at the PC surface and extending to 
O.SRic. The pair cascade simulation code that we use is based on one developed and described 
in several previous papers (Daugherty & Harding 1996, Harding et al. 1997, MHOS). The 
code follows a primary electron along a number of magnetic field lines above the PC, starting 
at rest and tracking its Lorentz factor from a given accelerating E\\ and energy losses due 
to CR and inverse Compton scattering. The paths of photons from curvature and inverse 
Compton radiation are traced to determine if they pair produce (by the one-photon process) 
or escape, at which point they are accumulated in a 3-dimensional array of energy and angle. 
Photon splitting, as treated in Harding et al. (1997) has not been included in the present 
calculation, and indeed is not important for most pulsars in the SG at high altitudes where the 
magnetic field is well below the critical value. Synchrotron radiation from several generations 
of secondary pairs is also traced, until all radiated photons escape the magnetosphere. As 
in MHOS, we adapt the primary acceleration part of the code to include only acceleration 
in the SG; thus we inject primary electrons only along the last open field lines at equally 
spaced azimuth about the magnetic pole. We have assumed a static vacuum dipole model 
for the magnetic field and thus do not consider the distortions to the PC footpoints or the 
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bending-back of field lines near the light cylinder that appear in retarded vacuum dipole 
models (e.g. Dyks et al. 2003). Since the pair cascades occur at relatively high altitude, 
above the pair formation front forming the inner bounary of the SG, where the E\\ is screened, 
the secondary pairs are not assumed to accelerate. With the E\\ of eq. (59), only primary 
electrons continue accelerating to high altitude within the SG. As shown in Section 3.1, the 
accelerating primary Lorentz factors become radiation-reaction limited by curvature losses 
at altitudes within and above the pair cascade zone (1.5 — 3) R. Thus, these particles will 
continue radiating well above the cascade zone with Lorentz factor given approximately by 
eq. (56). 

Figure 1 shows an observer angle-phase {( — (j)) plot of escaping photons above 100 MeV 
for pulsar inchnation angle x = 45° and parameters P — 0.033 s and Bq — 8x 10^^ G of the 
Crab pulsar. One can see the pattern of field lines due to the discrete spacing of azimuthal 
injection points of the primary electrons. The intensity of radiation along each field line is 
shown as a grey scale. One of the important features to note is the ring of enhanced emission 
around each magnetic pole, which is the hollow emission cone formed by the pair cascade 
radiation. This component terminates at altitudes around (3 — 4) R. Another important 
feature are the lines of caustic emission trailing out beyond the pair cascade rings from 
each pole. The caustic emission comes from electrons radiating curvature radiation between 
altitudes of (0.1 — 0.7) Ric along the trailing field fines. As first noted by Morini (1983), who 
considered only the caustic emission from one pole, and recently by Dyks & Rudak (2003), 
who considered emission from both poles, the caustics form because the positive phase delays 
from different altitudes along the traifing dipole field lines are nearly completely cancelled by 
negative phase shifts caused by aberration and time-of-flight. Since our code includes such 
relativistic effects, the two-pole caustic emission seen by Dyks & Rudak (2003) is a feature 
of extended SG acceleration. 

Pulse profiles result from slicing the phase plot of Figure 1 at constant observer angle, 
(. It is apparent that observers having 30°<C~70° and 120°<C~160° will see pair cascade 
emission, in which case either a single peak or two peaks with separation dependent on 
observer angle and inclination result. Observers having 60°<C^120° will cross lines of 
caustic emission from both poles, resulting in two narrow peaks in the pulse profile. Since 
the caustic emission lines are slanted with respect to each other in the phase plot, the peak 
separation will vary with observer angle and will, in general, not be 180°. Figure 2 shows 
a profile for C = 113°, in which there are two peaks separated by about 140° with a higher 
level of emission between the peaks, about what is measured for the separation of the peaks 
in the Crab and Vela high-energy pulse profiles. Note that the closest approach to one of 
the magnetic poles occurs at phase —180°, where an observer could also see the edge of a 
radio conal emission component. This could explain the position of the precursor in the Crab 
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pulsar, or the phase of the radio pulse in the Vela pulsar, both of which lead the double 7-ray 
peaks. An identical profile would be seen at observer angle 67°, which is approximately the 
viewing angle inferred from the Chandra image of the inner Crab nebula in X-rays (Weisskopf 
et al. 2000). The caustic emission occurs on the unfavorably curved field lines, i.e. those 
curving away from the rotation axis at observer angles X < C < 180° ~ X- The inclusion 
of the frame- dragging component of the potential, which is the only one that survives at 
high altitude and is not dependent on magnetic azimuth, is thus critical to allowing caustic 
emission. At certain (very limited) observer angles, it is possible to view both cascade and 
caustic emission, in which case there would be three or four peaks in the profile of varying 
relative intensity. 

Acceleration and cascade emission along field lines interior to the SG ('core' emission) 
is also expected, but has not been included here. Since the pair cascades on the interior 
field lines occur much closer to the NS surface, the emission would form a hollow cones of 
smaller radius in the phase plot. For observer angles crossing near the magnetic pole, such 
core emission would produce flux in the interpcak region of the profile. For observer angles 
crossing the caustics, the core emission would not produce much additional flux in the proflle, 
at least at high energies, since the electric fleld is screened at low altitude on the interior 
fleld lines and primaries do not continue to accelerate above the cascade region as they do 
in the SG. 

In the present simulation, we have included only CR of the primary electrons accelerat- 
ing at high altitudes, which forms the caustic emission. As noted earlier, they have reached 
radiation- reaction limit so that their CR spectrum will be quite hard (photon index —2/3). 
This spectrum is much harder than the observed spectra of 7-ray pulsars. However, the 
primary electrons will also be scattering sources of soft photons, such as radio, optical and 
infrared. We are investigating the possibility that inverse-Compton scattering of radio emis- 
sion along the same fleld lines, forming the conal component in phase with the high-energy 
peaks, may account for the spectrum of the Crab pulsar. 

4. Discussion and Conclusions 

In this paper we extend our previous study (MH03) of pulsar SGs to investigate the 
regime of acceleration of primary electrons and high-energy emission to very high altitudes. 
Incorporating the effect of cross-fleld motion of charges near the light cylinder on the dis- 
tribution of accelerating potential within the entire SG, we derived the explicit expressions 
for the accelerating electric fleld in the space-charge-limited flow approximation. We have 
modeled the particle acceleration up to high altitudes and generation of high-energy emission 
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within the SG, illustrating the resulting high-energy radiation pattern and light curves for 
the case of the Crab pulsar. The primary (space-charge-limited) flow within the SG becomes 
radiation-reaction limited at high altitudes, such that the energy gain from acceleration is 
nearly compensated by the CR losses. The resulting emission pattern exhibits both the hol- 
low cones centered on each magnetic pole (the corresponding phcnomenological model was 
first discussed in connection with radio emission by Komesaroff 1970, Radhakrishnan 1969, 
and Radhakrishnan & Cooke 1969), due to pair cascades on the inner edge of the SG, as 
well as the caustics along the trailing field lines at high altitude noted in previous studies 
(Morini 1983, Romani & Yadigaroglu 1995, Cheng et al. 2000, Dyks & Rudak 2003). The 
extended SG acceleration allows the caustic emission to be viewed from both magnetic poles 
simultaneously, naturally producing the widely-spaced double-peaked profiles seen in many 
high-energy pulsars. 

Note that in a vacuum NS magnetosphere with the strong magnetic field prohibiting 
transfield motion, the GJ space charge may serve as a source of electric field accelerating test 
particles (e.g. electrons/positrons) along the magnetic field. In the SG geometry with space- 
charge-limited flow (nearly GJ near the surface), the occurrence of a steady-state cross-field 
drifting at very high altitudes tends to reduce the GJ space charge (see condition [15] and 
related discussion) within the SG, maintaining the SG boundaries as equipotential surfaces. 
The reduction of the GJ space charge is accompanied by a steady-state distribution of the 
effective flux, AF, within the SG. Our model requires that the source of a steady-state 
accelerating electric fleld within the SG establishes at relatively low altitudes (e.g. over a 
few stellar radii above the PC) and remains constant along the magnetic stream lines up to 
very high altitudes. 

Whether or not a steady-state particle flow can be achieved on a given fleld line will 
depend on the relative importance of the terms oc cos x and oc sin x cos (ppc near the NS 
surface, which depends on the pulsar inclination angle x, the magnetic azimuth cos0pc a-nd 
the PC opening angle Oq- If the term oc cosx dominates, then E\\ accelerates the same sign 
of charge at all altitudes and steady-state flow can be maintained. However, if the term 
oc sin X cos begins to dominate near the surface for large x on unfavorably curved fleld 
lines (cos0pc < 0), then E^^ changes sign and may prevent a steady-state flow of charge (or 
any charge flow) along that fleld fine (see Section 3.2.1). 

Even if E\\ does not change sign above the PC, the possibility we have explored in this 
paper of a steady-state regime of acceleration of primaries from the PC surface up to the 
light cylinder implicitly assumes a free supply of charge from the NS. We do not exclude 
a scenario with non-stationary or insufficient supply of primaries from the PC surface, in 
which case there could be a regime of intermittent charge flow. Note that we have briefly 
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discussed a similar possibility in the case of acceleration of primaries over the main part of 
the PC (see HMOl). 

Outer gap models (e.g. Cheng et al. 1986, Romani 1996), that assume vacuum gaps 
form along field lines above the null charge surface, require little or no current flow into 
the gaps from outside. In the model presented in this paper, steady current flow can exist 
from the PC to the light cylinder along even some unfavorably curved field lines where outer 
gaps are assumed to exist. In this case the space-charge limited current flow of nearly GJ 
surface value, combined with cross-field particle motion, would prevent the formation of an 
outer gap. In the case of very high inclination angle, where cither non-steady current flow 
or insuflicient charge supply from the surface exists, outer gaps may still form (see Section 
3.2.1). In this case, however, the radiation which forms caustic emission above the null 
surface can be viewed from only one pole (Romani & Yadigaroglu 1995, Cheng et al. 2000). 
The geometry of the extended SG emission is thus clearly distinct from that of the outer gap 
emission and will have observable consequences, as has been discussed by Dyks & Rudak 
(2003) and Dyks et al. (2003) and will be detailed in future pubhcations. 

We would like to thank Vladimir Usov, Kouichi Hirotani, Bing Zhang and Jarek Dyks 
for valuable discussions and comments on the manuscipt. We also acknowledge support from 
the NASA Astrophysics Theory Program. 
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Fig. 1. — Map of emission above 100 MeV as a function of phase and observer angle ( 
with respect to the rotation axis for a pulsar with magnetic inclination angle of % = 45°, 
period P = 33 ms and surface magnetic field Bq — S x 10^^ G. Magnetic poles are located 
at (0 = 0°, C = 45°) and (0 = 180°, C = 135°). 
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Fig. 2. — Pulse profile for emission above 100 MeV for the case of Figure 1 and observer 
angle C = 113°. 



